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Abstract 

The Landau distribution as well as its first and second momenta are well suited for describing 
the energy loss of charged particles traversing a thin layer of matter. At present, just rational 
approximations and asymptotic expressions for these functions were obtained. In this paper we 
present a direct calculation of the integral representation of these functions obtaining perturbative 
and nonperturvative solutions expressed in terms of fast convergent series. We also provide a simple 
numerical algorithm which allows to control speed and precision of the results. The testing runs 
have provided, in reasonable computing times, correct results up to 13-14 significant digits on the 
density and distribution functions and 9-10 on the first and second momenta. If necessary, this 
accuracy could be improved by adding more coefficients to the algorithm. 
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I. INTRODUCTION 

It is well known that the Bethe-Bloch formula[lj] describes the average energy loss of 
charged particles when travelling through matter, while the fluctuations of energy loss by 
ionization of a charged particle in a thin layer of matter was theoretically described by 
Landau[|. This description ends with a universal asymmetric probability density function 
characterized by a narrow peak with a long tail for positive values. This tail towards positive 
values comes from the small number of individual collisions, each with a small probability 
of transferring comparatively large amounts of energy. An integral representation of the 
Landau probability density function reads 

1 f°° 

(p(x) = — y y sin (ny) exp (— xy) dy (1) 
7T Jo 

with x > 0. 

Other functions such as the corresponding distribution function $ (0, x) and the first and 
second momenta $(l,x) and $(2,x) of the density function truncated on the right-hand 
tail defined through the general formula 



$ (n, x) 



(n-l)(n-2) , "(3-w) $ / 0ja .x 



c n <p (c) dc (2) 



oo 



2 2 

are also needed in order to fit a truncated Landau distribution $ (0, x) to measurements or 
to simulated energy-loss data 

It is clear that a direct calculation, by numerical integration, of expressions Q]) and 
(J2j) besides of being of a very slow convergence could be used only to provide values for a 
subsequent construction of approximate formula and for testing purposes. 

In the seventies some programs for the numerical computation of the Landau distribution 
were presented J|. There, the main emphasis was put more on the speed of the calculation 
algorithms than in the precision obtained, that was in general not good enough. In the 
eighties a new package of programs^] arouse. These programs improved the precision of the 
values obtained for JH) and (J2J) by using piecewise functions as rational approximations and 
considering their asymptotic behavior. One should mention also Ref. 6 were more accurate 
calculations of the Landau density were presented. For more recent articles see Ref. 7. In 
view of the still high computational time required, the precision for the distribution function 
was only up to 6-7 significant digits. At that time, this was thought to be enough for the 



applications in experimental physics but their use in other fields of application remained 
excluded. 

For practical uses and under the standard prejudice that equation {TJ cannot be rewritten 
into an analytic form, the following approximate expression was suggested^] 



This analytic approximation, known as the Moyal function, is in fact an extreme value 
distribution 2|. It provides values that are not entirely coincident with those coming from 
the original expression (|Tj). in particular in connection with its maximum value and its 
asymptotic behavior. Nevertheless, being quite simple, it is widely used in different fields. 
See Ref. 10 for examples. 

To overcome the aforementioned difficulties, new and refined approximations, combin- 
ing speed and accuracy in the evaluation of these functions, are required to improve the 
agreement between measured and calculated energy loss distributions. 

Under this motivation, we have carried out a direct calculation of the integral represen- 
tation of the Landau density, the distribution function and the two first momenta of the 
density. We have obtained perturbative and nonperturvative solutions expressed in terms of 
fast convergent series. To achieve this, we have used analytical continuation properties and 
the Laplace transform of the Landau density. We have also developed a simple numerical 
algorithm which allows one to control speed and precision of the results. 

In what follows we summarize the main steps for obtaining the explicit analytic expres- 
sions for the Landau probability density JTJ and the momenta (J2J) . An appendix with details 
of the calculation is annexed. In addition, a brief description of the algorithm which eval- 
uates numerically these functions is presented. Finally, a numerical test together with a 
critical analysis of the results is provided. 



4> (x) ~ exp — - (x + exp (— x)) 




(3) 
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II. NONPERTURBATIVE SOLUTIONS 



A. Landau density 

For computational purposes, it is of interest to consider the following auxiliary integral 

1 f°° 

Q(x,b) = - y y exp[-y(x + b)}dy. (4) 
Jo 

Clearly, once its solution is known, it is necessary to perform the analytic continuation 



(x) = Im [Q (x, b = —in)] 



(5) 



to obtain the Landau density. 

By using standard representations of the functions in the integrand of Q(c,b), one can 
write 



oo 

Q(x,b) = -J2 



) k d k 



TT * — ' k\ da k 

k=0 



T(k + a + 1) 
(x + b) 



k+a+1 



a=0 



where T is the gamma function(ll|. Using now the identity 

^k k 

- [F (a) G (a)] = J 

r=0 \ r 



[F (a) G (a)] = £ [ ] ^-[F (a)] — [G (a)] 



da k 



da k ~ 



da r 



and defining the coefficients 



.4 



kr 







k—r 



da k ~ 



T (k + a+1) 



(6) 



a=0 



one obtains 



fc=0 ' r=0 \ r / 









[r. 


1 


{x + b) k+1 _ 



[In (x + b)} 1 



Finally, performing the analytic continuation (JHJ), one obtains the exact analytic expression 
for the Landau density 



[x 



n oo / xfc k 

ii ^— ' k\ ^— ' 

fc=0 r=0 



r 



) r A kr lm 



[In (x — ivr)] r 

"7 • \k+l 

(X — 171) 



(7) 
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Notice that this expression can be analytically continued to x < values. However, we 
obtain below a formal solution also for this case, since the numerical evaluation of eq.(0) for 
negative values of x makes the series divergent. 

B. Landau distribution and momenta 

In order to compute the functions defined in Eq.Q we use the Laplace transform of the 
Landau density, namely 

exp (— pc) (f> (c) dc = p p ; p > 0. 

3 

The derivative of this relation with respect to p, gives rise to a new and convenient 
expression for computing the Landau distribution and the two first momenta of the density 
as follows 



$ (n, x 
lim 

p -> 0+ 



ii 



( B -l)(n-3/2) . 

P 



(8) 



c n e- pc 6(c) dc 



Replacing now the Landau density (JZJ), into this last expression, we see that the calculation 
is reduced to the computation of 



k=0 r=0 y T J 



kr 



lim 

p -> 0^ 



Im 



+oo 



n .- pc [ ln ( C -^)] 



c e 



(c — in) 



fc+i 



dc 



that means the integrals 



(9) 



lim 

Bt (P, a;) = ^ Im 

p -> 0+ 



c n e - pe [ln(c-^)r dc 



(c — m) 



fc+1 



(10) 



Once these integrals are evaluated, an analysis of their eventual divergences and the 
necessary cancellations are in order. 
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1. Landau distribution 



Let us solve first the case n — 0, the Landau distribution function. Notice that in this 
case all the integrals in expression ([TO are well defined and consequently the interchange 
between the limit operation and the integral is licit. In this way we obtain 



B° 00 (0,x) = lm 



1 



C — ITT 



-de 



7i [1/2 — arctan (x/tt) /tt] 



and 



Bl (0, x) 



lim 

b — > OCT 



Im 



ln(6— in) 



u r exp (—ku) du 



\n(x— I7r) 



(11) 



(12) 



A further change of variables to v = uk, reduces the integral to an incomplete Gamma 
function [ijj. Finally, after taking the limit, one obtains 



B° kr (0,x) 

Going back to equation (JHJ) one gets 



Im [r (r + 1, k\n (x — in))} 



11 ^OO/xfcfc 

G (0, x, 0) = - - - arctan (z) + ~Y1 j^F 



2 7T 

so that the Landau distribution results 



7T/ 7T •* — ' k\ 

k=l r=0 \ r 



A kr B° kr (0,x) 




11 (X 
— I — arctan — 
2 tt Vtt 



k r+1 



Im [r (r + 1, A; In (x — i7r))] 



(13) 



£. Momenta 

Going now to the analysis of the momenta, we found that in the general expression 
there appear divergences that finally cancel out, but that have to be treated with care. We 
postpone to an Appendix some details of the lengthy calculation involved and present here 
the final results, both for the first and the second momenta. They read 



$ (1, x) $ (0, x) = -1 + 7 + - In (x 2 + tt 2 ) 
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and 



I 1 

+- V (-) r A lr [nRe [T (r + 1, In (x - in))] + 

7T Z — ' 

r=0 

1 1 1 1 1 

[l — r + r In (x 2 + 7r 2 )] [?t/2 — arctan (x/V)] Im (271") 1 



m=0 



E E ( ^ ) n7 ] \ A X [r (r + 1, (fc - m) In (x - ztt))] 



fc! 

fe=2 r=0 



r / (k — m) 



(14) 



7T 



$ (2, x) $ (0, x) = 1 + x + - [7r - 2 arctan (x/7r) 



+ (l- 7 ) 



a; 



1 - 2 



1 

7T 

1 

7T 



X 2 + 7T 2 ) 

x 2 [x [arctan (x/tt) — f] + § In (x 2 + 7r 2 ) + 7r] 

X 2 + 7T 2 



7 + - In (x 2 + 7T 2 ) 



7T 



271 + 77T + 7 7T + — 

o 



7T 



[in (x 2 + 7T 2 )] 2 — 2x arctan (x/tt) — — + n arctan (x/tt) 



TT 



1 2 

-V 

2tt ^ 



^ i I 27T 



E 





(-) r A 2r Im 


( r J 







r+l 



' \ m / (2-m) r+ 

1 2 
— Im > 

TT 



[In (x — 271")] 

r + l 



r [r (r + l,(2-m) In (x - in))] 



2 1 / A 2-m (~) k 

m) { ) ^ 

m=0 \ / k=3 



E 



.4 



At 



\r+l 



[r (r + 1, (A; — m) In (x — wr))] 



(15) 



r=o \ r I (k-m) 

being 7 the Euler's constant [Uj. Notice that the limit x — > 00 indicates that these distri- 
butions, even if normalized to 1, do not have well defined momenta. 
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III. PERTURB ATIVE SOLUTIONS 



A. Landau density 

For the sole purpose of presenting solutions in an as simple as possible way, we define 
two auxiliary functions Mj (x) and (x) . First 



M j (x) 



D 



jnkp 



7T L — ' n\ 

n=0 k=l 



p=Q 



(2p + 1)!4p' 



(16) 



where the coefficients Dj n k p are given by 



D 



jnkp 



d 2 P 
Q u 2p 



(u + k) 



(—j—l—u—k) n • I ^ 

v J 'a sin to, 



(17) 



«=l/2 



Note that the expression (fTo|) is the pertubative solution of the integral 



1 f°° 

Mj (x) = — / y-v-i- sin (vry) exp (—xy) dy. 
* Ji 

This solution is simply obtained by dividing the integration range in consecutive unit 
segments and expanding afterwards the exponential and trigonometric functions in the in- 
tegrand. 

The other auxiliary function is given by 



Tiji (x) = 



— E 



7r L — ' n\ 

n=0 



'tijln 



(18) 



being the coefficients Uji n given by the following expression 



tijln -2^ (2m + l)\ 



m=0 



E 



(k+j)\ 



1(1-1) /2 
2m + 3 J f^ k\(k + n + i + 2m + l + l 



The expression |THj) is the perturbative solution of the integrals 



7T 



y cos (7T7/) exp (—xy) In 7 (y) dy, — / y y+l sin (717/) exp (—xy) ln J (y) 



7T 



and 



- / 2/ 

7T 



COS (7T7/) 



sin (%y) 



ivy 



exp (—xy) In 3 (y) dy 
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for 1 — 0,1 and 2 respectively, which is obtained by expanding the exponential 
exp (— yln (y)) , and the other functions appearing in the integrand, in powers of their argu- 
ments. 

As a consequence, the Landau density is obtained by dividing the range of integration in 
expression ((TJ) into [0,1] x [1, oo]. Thus, we get 

<\> (x) = Tooi (x) + M_! (x) . (19) 



B. Landau distribution and momenta 

The computation of these functions follows the steps done for the case of the nonpertur- 
bative regime. In this case, the integrals appearing in Eq.fjHJ) are now computed from the 
following relationship 



/ c n e^ c (c) dc=-^2\ ./!•'•" J / ~, tt+t sin (Try) exp (-xy) dy. (20) 

•/* 71 j^o \ j J Jo (p + v) 

In the limit p going to zero, Eq.JBJ evaluated at n = becomes 

1 r 00 y-y 

$ (0, x) = 1 / - — ■ — r sin (Try) exp (-xy) dy. (21) 

n Jo (p + y) 

Let us observe that this integral is well defined allowing us to evaluate it at p — 0. In 
this way, we obtain the following expression for the Landau Distribution 

1 f°° 

$(0, x)—l / y~ v ~ x sin (rry) exp (-xy) dy. (22) 

71" Jo 

Analogously, we get for n = 1 and n = 2 

$ (l,a;) $ (0,x) = -1 - In (p) / y~ v ~ x sin (7n/)exp (-xy) dy (23) 

Jo 

1 r°° y-y 
/ ~ 72 sin (iry) exp (-xy)dy 

71 Jo (p + y) 

and 
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$ (2, x) $ (0, x) = 1 + 3 In (p) + - + In 2 (p) 



/ -3 sin (Try) T (3, z(y + p))rfy (24) 

(p + y 



_i r 

^ Vo (p + y) £ 

corresponding to the first and second momentum respectively. 
We can conveniently rewrite these last two equations as 



1 f°° y~ y 

$ (1, x) $ (0, x) = -1 - In (p) - x [1 - $ (0, x)] / — g sin (7ry) exp (—xy) dy (25) 

^ (p + y) 



and 

1 



$ (2, x) $ (0, x) = 1 + (3 + 2x) In (p) + - + In 2 (p) + 2x + x 2 [1 - $ (0, x)} 

V 

2 Z" 00 y-v 

+2x$ (l,x) $ (0,x) / 3 sin (7ry) exp (-x (y + p))dy. (26) 

71 Jo (p + y) 

By simplicity we also divide the range of integration of the integrals present in (|22|) . (|2o|) 
and (|26|) as [0, 1] x [1, oo] such that we can easily evaluate them in terms of the two auxiliary 
functions previously defined. In particular, all the integrals defined in [1, oo] are well defined 
and consequently they can be evaluated at p = before being computed. The divergent 
behaviour in the parameter p of the integrals defined in [0, 1] is extracted by integrating 
them by parts. Then, the remaining finite parts are also expressed in terms of the auxiliary 
functions. 

In the case n = the integral defined in [0, 1] is well defined and we do not need to 
integrate it by parts. Consequently, Eq. (|2*2T) provides directly the final expression for the 
Landau Distribution 

$(0,x) = l-Mo(aO-r_ 10 i(aO (27) 
which for the first and second momenta one has 



= _x[l_-$(0,x)]+M 1 (z) 



$(0,x) 

T-121 (x) + (l + x) T_ m (x) + ttT_ii 2 (x) 
$(0,a;) 



(28) 
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and 



$(2,x) 



1 + 3x + x 2 [1 - $ (0, x)} + 2x$ (1, x) $ (0, x) 
$(0,x) 

exp (-x) [1 - tt 2 /3] - 2M 2 (x) - / (x) 



$(0,x) 

respectively, where / (x) is given by the expression 



(29) 



Kx) 



-if 

-7rT_ 



-7T T ( 



7T 2 

T 02 i (x) + 2 (1 + x) Ton (x) - y^ooi (s) 
-ii2 (x) + (1 + x) T_m (x) - 7rT_i i (x) 



+ 



vr 2 

030 0) + 3 (1 + x) T 020 (x) - — (l+x) T 000 (x) 



2(l + x) 2 - y 



^oio (x) 



(30) 



IV. ALGORITHM 

The nonperturbative expressions 0, (fl3|) . |H| and (fl~5f are used for x > 3 , and the cor- 
responding perturbative expressions (IT9l) . (|27|) . (|28| and (|29| for |x| < 3. These expressions 
are the kernel of our algorithm. 

In order to compute the nonperturbative expressions we have provided the algorithm 
with the following numerical coefficients 



\ k+r 



A 



Bkr ~ r! (Jfc-r)! (31) 
where A kr is given by (|EJ). These coefficients were computed for k > r by the following 
relationship 



.4 



kr 




r 
P 



daP 



r(k + a + l)V(k-r-l-p,k + a + l) 



(32) 



a=0 



^ (n, z) being the polygamma function|ll| and A kk = T {k + 1). 

We computed these coefficients through the following numerical recursive algorithm 



11 



forj = 0,(k-r)^A(k,j) = Yl { 3 1 ] A(k,p)V (j - 1 - p,k + 1) 

p=o\P J 

where A(k,0) = F (k + 1). Consequently A kr = A(k,k — r). Replacing these coefficients 
into Eq. pTf we get the numbers B kr . This algorithm was worked with MAPLE software. 

Taking into account that our results provide excellent approximations, considering only 
a few terms of the corresponding series, we have just calculated the first 5050 coefficients 
corresponding to values of k = • -99 and r — • -k 

The computation for the functions $ (0, x), <E> (1, x) and $ (2, x) involves, in addition, the 
evaluation of the Incomplete Gamma function with complex argument and integer order 
T (n, z). It is accomplished by using the following expression [n| 

n-l „ 



T (n, z) — (n — 1)! exp (— z) — - n > 1. 



On the other hand, the perturbative expressions are computed by providing the algorithm 
the following numerical coefficients 

D2jnkp = (2p+"l)!4P (33) 
where Dj n k p is given by the Eq. (jl7|) . These coefficients were numerically computed for 
< j < 3, < n < 26, 1 < k < 20 and < p < 14 using MAPLE software. 

In addition we have defined auxiliary functions to compute the expressions JUJ), (fT8|) and 
JUJ. 

The data files containing the aforementioned numerical coefficients (pTTJ) and fHHl) can be 
found in the following link: www.fisica.unlp.edu.ar/download/landau. 



V. NUMERICAL RESULTS 

The test of our algorithm is based on the numerical integration of the expressions (JTJ) 
and (121) using MAPLE software. We found that the algorithm allows to obtain results 
with a precision up to thirteen to fourteen significant digits for the Landau density and 
distribution and nine to ten for the first and second momenta. If necessary, higher accuracy 
can be obtained by adding more coefficients B kr and D2j nkp into the algorithm. 
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The set of points chosen to evaluate the aforementioned functions includes those values 
presented in Ref. 5 where the precision achieved certainly depends on the convergence of 
the expressions obtained. The corresponding numerical solutions for these functions are 
tabulated in the Tables 1,2,3 and 4. We also included in these Tables an estimation on the 
order of the series associated with our solutions by which a lower precision in the numerical 
results is obtained. Specifically, we have computed |A/(x)| = \f s+ i(x) — f s (x)\ < 10~ d , 
where / (x) represents any of the following four functions; 4> (x), $ (0, x), $ (1, x) and $ (2, x), 
and "s" represents the number of terms needed by the sum over the dummy variable k (Eqs. 
©, fTCty ; (TT4T) and fT5|0 or n (Eqs. ljTB)) and (fll?jl ) to obtain a precision of "d" significant 
digits in the nonperturbative or perturbative regimes respectively. 

Indeed, this information provides insight on how the convergence of our solutions depend 
on the point where the function is evaluated and on the required precision. In particular, 
we note that faster evaluations correspond to points localized closer to the origin or to the 
positive tails. 

For getting some computing time estimation, we developed a partially optimized program 
written in Fortran 77 which is based on the above described algorithm. We found that the 
typical running time to compute, with a given accuracy, the four functions over the complete 
set of points is, in general, of the order of the fraction of a second in a common PC processor. 
For instances, the total number of CPU seconds that the process spent in user mode is around 
0.205, 0.197 and 0.185 for d equal to 12, 9 and 6 respectively. 

VI. SUMMARY AND DISCUSSION 

The Landau distribution as well as its first and second momenta are used to describe the 
energy loss of charged particles by ionization, particularly within the area of the high en- 
ergy physics. At present, just rational approximations and asymptotic expressions for these 
functions were available. Numerical algorithms using these approximate solutions provide 
a fairly fast computation. In this paper we presented a direct calculation of the integral 
representation of these functions obtaining perturbative and nonperturvative solutions ex- 
pressed in terms of fast convergent series. These solutions were obtained by using analytical 
continuation properties and the expression of the Laplace transform of the Landau density. 
They have been shown to be very flexible, allowing to combine accuracy and speed to adapt 
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our approach to the computation of these functions to many fields of application. We re- 
mark that with only few terms of our series and a fast evaluation of them in all the range 
of its variables, we obtain excellent approximations. In fact, we found that our algorithm 
provides in reasonable computing times correct results up to 13-14 significant digits on the 
density and distribution functions and 9-10 on the first and second momenta. This is clearly 
a substantial numerical advantage over the existing programs mentioned above. If necessary, 
one can improve that precision by including further coefficients into the algorithm. 

Future work involves the obtention of approximate analytical expressions for the single 
collision cross section and the electro dynamic properties of the detector medium These 
quantities, which can be computed in terms of the Landau distribution, could preserve the 
main features of our solutions. In particular they could be similarly expressed in terms of 
fast convergent series. This will surely improve not only its numerical evaluation but it will 
also allows to produce suitable fits of the truncated Landau distribution to measured or 
simulated energy-loss data. In fact, one should obtain an approximate analytical expression 
for the dielectric function (or equivalently the oscillator strength function) which is well- 
defined over the complete range of energy. This should remove the presence of non-physical 
high energy behaviour introduced in Ref. 12 by using a different approximation for the 
Landau distribution. 
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Appendix A: TECHNICAL DETAILS 

We start by analyzing the expression (jTUj) B% r (p, x) where we cannot, a priori, interchange 
the operations of limit p — > + and the integration. In those cases we are interested in the 
obtention of their finite analytic contribution. Their divergent terms cancel with terms like 
pP [n (n — 1) (n — 3/2) /p + (1 + In (p)) n ] coming from Eq. (JHJ). Then we solve the other 
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integrals where the interchange of operations simplifies the calculation. 

In the case n = 1 the divergent integral is just the first one which is given for k 
r = 0, namely 



= and 



lim 

B oo (P, x)= Im 
p -> 0+ 



cexp (—pc) 



C — %1T) 



-de 



Integrating by parts and expanding for small p we obtain 



7 + - In (x 2 + vr 2 ) + In (p) 



+ C(p). 



(Al) 



The other integrals corresponding to (k > 1) can be rewritten, after a change of variables, 



B\ r (0, x) = Im I 

m=0 \ m 



1— m 



lim 

6 — > oo^ 



ln(6— iir) 



u r exp [ — u (k — m)} du 



\n(x—iir) 



Taking into account that these integrals can be expressed in term of the incomplete Gamma 
function, we obtain: 

• for k > 1 



m=0 \ m 



in) 1 -™ ^—^ [T (r + 1, (k - m) In (x - itt))] (A2) 



(k — m) 



for k = 1 



£j r (0, x) = nRe [r (r + 1, In (x - iir))} 



+ 



— [l — r + r In (x 2 + 7r 2 )] [7r/2 — arctan (x/tt)] 



(A3) 



Replacing Eqs JAljl . (|A2jl and (|A3jl into Eq.^ we check that the logarithmic divergent 
terms finally cancel and the remaining part of the integral provides the solution given by 
expression (jH|) . 

In the case n = 2 the divergent integrals are for k = and = 1 and they read 



#oo (P. 



lim 



Im 



-OO 1 

^2 



c exp (—pc) 



-de 



(c — in) 



and 
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2 ^ m 

p -> 0+ 



+00 



c 2 exp (— pc) 



[In (c — i7r)] T 



C — 27T 



respectively. Integrating by parts and expanding for small values of p we obtain 



for k = and r = 
for = 1 and r = 



1 7T 

x — — [tt — 2 arctan (x/71")] 



+ C(p) 



(A4) 



5 io (P> a;) = 7T 



(x 2 + 7T 2 ) 



1 - 2 



7 + In (p) + - In (x 2 + 7r 2 ) 



+ O (p) (A5) 



for k = 1 and r = 1 



£11 (P, 



[x [arctan (x/tt) - | ] + | In (x 2 + tt 2 ) + tt] 

X 2 + 7T 2 



— 27T + 77T + 7 7T 



TT ' 



- [In (x 2 + 7r 2 )] 2 -2x 



71 



arctan (x/ir) — — 



7l 



71 



arctan (x/71") — — 



1 2 



+7T (1 + 2 7 ) In (p) + vr [In (p)] 2 + (p) (A6) 
The others integrals (k > 2) can be solved in the same way. The integrals can be written as 



£ 2 r (0,x)=lm£( 2 W- Um 



that allows one to obtain 



ln(fe— itt) 



w r exp [ — u (k — m)} du 



ln(x— in) 



for k > 2 



5 2 r (0, x) = Im 



m=0 



* H 2 ~ m 7j ~T7+T P (r + !,(*- "0 h 0* " **))] (A7) 



for A; = 2 
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B 2 2r (0,x) = Im 




2-m 



T [[ln(x-^ 7 r)r +1 ] 



[r (r + 1, (2-m) In (x - in))} 



m=0 



(A8) 



(2 - m) r+ 

Replacing now Eqs.JSU), (EH), (jMJ, flSH) and l(X8]) into Eq.© we explicitly check that 
the divergent terms 1/p, 3 In (p), and In (p) 2 finally cancel and the rest of the integral provides 
the solution given by the expression for the second moment. 
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Table 1 

Numerical results for the Landau density. The integers given in columns 3-5 represent the 
number ofterms needed to achieve an accuracy of 12, 9 and 6 significant digits respectively. 



X 


<j>{x) 


|A0| < 10- 12 


|A0| < 10~ 9 


|A0| < 10- 6 


-3 


3.658377588806938£ - 002 


20 


16 


12 


-2.5 


9.637807619860883£ - 003 


18 


15 


11 


-2 


4.398547840053790£ - 002 


17 


14 


10 


-1.5 


0.100550751729020 


15 


12 


9 


0.4 


0.168703915408766 


10 


8 


6 


0.5 


0.165232275482657 


11 


8 


6 


1.5 


0.124221094071471 


16 


13 


9 


3 


7.424765459960213E - 002 


21 


18 


14 


4 


5.326865722328474E - 002 


29 


17 


8 


5 


3.916341958128983£ - 002 


21 


13 


7 


8 


1.809028394181073E - 002 


10 


9 


7 


10 


1.197648738878854E7 - 002 


11 


9 


6 


15 


5.401336928010399£ - 003 


10 


8 


6 


20 


3.004979395692253£ - 003 


10 


8 


6 


30 


1.298671406775338^ - 003 


9 


7 


5 


60 


3.082734332814474E - 004 


8 


6 


4 


100 


1.076112240078372E - 004 


7 


6 


4 


500 


4.085712883475095£ - 006 


5 


4 


2 


1000 


1.012057179716136^-006 


5 


3 


2 


10000 


1.001659320662497^ - 008 


3 


2 
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Table 2 

Numerical results for the Landau distribution. 



X 


$(0,:z) 


|A$| < 10~ 12 


|A$| < 10~ 9 


|A$| < 10~ 6 


-3 


-1.607879630341813E - 003 


20 


17 


13 


-2.5 


1.961867909885173^-003 


19 


16 


12 


-2 


1.409435793315372E - 002 


17 


14 


11 


-1.5 


4.982427497490971E - 002 


15 


13 


9 


0.4 


0.356581474467521 


10 


8 


6 


0.5 


0.373280562479680 


11 


9 


6 


1.5 


0.518343029194760 


16 


13 


10 


3 


0.664205997524218 


21 


18 


14 


4 


0.727271852143877 


26 


15 


8 


5 


0.773026779942834 


18 


12 


7 


8 


0.853459958739293 


12 


9 


7 


10 


0.882938865913564 


11 


9 


7 


15 


0.923293944585224 


11 


9 


6 


20 


0.943462646851257 


10 


8 


6 


30 


0.963250172910313 


10 


8 


6 


60 


0.982280577955328 


9 


7 


5 


100 


0.989571029163290 


8 


6 


5 


500 


0.997976630597665 


6 


5 


3 


1000 


0.998993481974151 


6 


4 


3 


10000 


0.999899912058608 


4 


3 


1 
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Table 3 

Numerical results for the first momenta of the density. 



X 


$(l,x) 


|A$(l,x)| < 10~ 9 


|A$(l,x)| < 10~ 6 


-3 


-2.94418823883569 


21 


18 


-2.5 


-2.67471836189293 


19 


16 


-2 


-2.25828545384356 


17 


13 


-1.5 


-1.86990165100088 


14 


11 


0.4 


-0.679860050176673 


9 


6 


0.5 


-0.629322356963745 


9 


7 


1.5 


-0.180019038237684 


13 


10 


3 


0.339504262821524 


18 


14 


4 


0.611170716736155 


15 


7 


5 


0.839831585464024 


12 


7 


8 


1.35514739591618 


9 


7 


10 


1.60809603494204 


9 


7 


15 


2.06972587987781 


9 


7 


20 


2.39437506213547 


8 


6 


30 


2.84447000017431 


8 


6 


60 


3.59268530830209 


7 


5 


100 


4.12998818980124 


7 


5 


500 


5.77818474253184 


5 


4 


1000 


6.47746787601202 


5 


3 


10000 


8.78657692470135 


4 


3 
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Table 4 

Numerical results for the second momenta of the density. 



X 


$(2,sc) 


|A$(2,x)| < 10~ 9 


|A$(2,x)| < 10~ 6 


-3 


8.66852146374132 


24 


21 


-2.5 


7.17765680763846 


22 


19 


-2 


5.14686389834850 


18 


15 


-1.5 


3.58354569903627 


16 


13 


0.4 


0.941718427356982 


9 


7 


0.5 


0.908679035129438 


10 


7 


1.5 


0.944186293458727 


15 


12 


3 


1.82656174148077 


21 


17 


4 


2.72090902630219 


18 


9 


5 


3.74974971267330 


13 


8 


8 


7.21517917106136 


10 


8 


10 


9.64861303605968 


10 


8 


15 


15.7892588337345 


10 


7 


20 


21.8615020920748 


9 


7 


30 


33.7144809006428 


9 


7 


60 


67.7261171683676 


8 


7 


100 


111.315504919124 


8 


6 


500 


525.970930432499 


7 


5 


1000 


1033.83849293122 


6 


5 


10000 


10066.7727522736 


5 


4 
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